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We address the viability of exact Yukawa unification in the context of general SUSY GUTs with 
universal soft-breaking sfermion and gaugino mass terms at the GUT scale. We find that this 
possibility is challenged, unless the squark spectrum is pushed well above the limits allowed by 
naturalness. This conclusion is assessed through a global fit using electroweak observables and 
quark flavour-changing neutral current (FCNC) processes. The problem is mostly the impossibility 
of accommodating simultaneously the bottom mass and the BR(_B — > Xs'y), after the stringent CDF 
upper bound on the decay Bs — > fJ-^ lJ-~ is taken into account, and under the basic assumption that 
the b ^ amplitude have like sign with respect to the Standard Model one, as indicated by the 
B ^ Xs£+£- data. 

With the same strategy, we also consider the possibility of relaxing Yukawa unification to 6 — r 
Yukawa unification. We find that with small departures from the condition tan/3 ~ 50, holding 
when Yukawa unification is exact, the mentioned tension is substantially relieved. We emphasize 
that in the region where fits are successful the lightest part of the SUSY spectrum is basically fixed 
by the requirements of b — r unification and the applied FCNC constraints. As such, it is easily 
falsifiable once the LHC turns on. 



I. INTRODUCTION 

The hypothesis of grand unification is able to ad- 
dress many of the unanswered questions of the Stan- 
dard Model (SM), like charge quantization or the 
quantum number assignments of the SM fermions. 
Moreover, augmenting grand unified theories (GUTs) 
with supersymmetry (SUSY) not only stabilizes the 
large mass hierarchy between the electroweak (EW) 
and the GUT scale but also leads to the possibility of 
exact gauge coupling unification. 

Although this unification works remarkably well 
in the minimal supersymmetric Standard Model 
(MSSM), in order to test the idea of grand unifica- 
tion one needs other, independent observables. A first 
candidate in this respect would be proton decay. How- 
ever, the absence of a signal at proton decay experi- 
ments constrains mostly non-supersymmetric GUTs, 
whereas the strong model-dependence of the proton 
decay rate contributions from dimension-five opera- 
tors makes it difficult to draw general conclusions 
[U, i, i, i, S, H, 0, i] on the viability of SUSY GUTs. 

A different way to test SUSY GUTs is by their 
predictions for the masses and mixings of the SM 
fermions. In this respect, SO(IO) is especially at- 
tractive because it unifies all quarks and leptons of 
one generation into a 16 representation of the gauge 
group, leading to the opportunity of a unified Yukawa 
coupling for the fermions of that generation. Though 
this is not phenomenologically viable for the two light 
generations, unification of the top, bottom, tau, and 
tau-neutrino Yukawa couplings might be possible if 
tan (3, the ratio of the vacuum expectation values of 
the two Higgs doublets, is close to 50. 



It must be taken into account, however, that the 
success of this Yukawa unification sensitively depends 
not only on tan/3, but also on the SUSY spectrum 
and parameters, because the Yukawa couplings are 
much more sensitive to weak scale threshold correc- 
tions than are the gauge couplings. In the absence 
of a clear signal in favour of supersymmetry or an a 
priori knowledge of the SUSY spectrum, one would 
need a set of additional observables to sufficiently 
constrain the allowed ranges for the SUSY spectrum 
itself, in order to test GUT predictions for fermion 
masses. It turns out that flavour-changing neutral 
current (FCNC) processes - loop-suppressed observ- 
ables that are highly sensitive to SUSY particle con- 
tributions - are especially suited for this purpose. 

In 0, an SO(IO) SUSY GUT model proposed by 
Dermi'sek and Raby (DR) fld\ , and featuring Yukawa 
unification, has been reconsidered in a global analysis 
in light of all the most precise data on FCNCs in the 
quark sector. While the model successfully describes 
EW observables as well as quark and lepton masses 
and mixings [13, [Hi , in Q it was found that the si- 
multaneous description of these observables and all 
the FCNC processes considered is impossible unless 
the squark masses are pushed well above the limits 
allowed by naturalness and within reach of the Large 
Hadron Collider (LHC). 

The aim of this Letter is twofold. First, we show 
that the problem pointed out in and mentioned 
above is a general feature of SUSY GUT models with 
Yukawa unification and universal sfermion and gaug- 
ino mass terms at the GUT scale, thus challenging 
the viability of these hypotheses, when considered to- 
gether. Our second aim is then to explore a possible 
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remedy, namely relaxing the hypothesis of Yukawa 
unification in favour of the less restrictive t — v and 
b — T Yukawa unifications. The departure from ex- 
act Yukawa unification can be quantified by the para- 
metric departure from the condition tan/3 ~ 50. As 
clarified below, this case will be relevant not only for 
SU(5), but also for SO(IO). This study will allow 
us to address the question whether a range of large 
tan/3 < 50 exists, where a successful prediction for 
the bottom mass and full compatibility with quark 
FCNCs are possible at the same time. 

II. YUKAWA UNIFICATION AND FCNCS 

It is well-known [H, [l^] that, under the assump- 
tions of a universal sfermion mass mie and a univer- 
sal gaugino mass mi/2 at the GUT scale, and with a 
positive /i parameter, Yukawa unification prefers the 
region in MSSM parameter space characterized by the 
relations 

-^o«2toi6, /i, mi/2<TOi6, (1) 

because they ensure a cancellation of potentially large 
tan /3-enhanced SUSY threshold corrections to the 
bottom quark mass which could otherwise spoil 
the Yukawa unification. Through renormalization 
group effects, these relations lead to an inverted scalar 
mass hierarchy (ISMH) [l^, i.e., light third genera- 
tion and heavy first and second generation sfermions. 

Relations ([T]), together with the large value of 
tan /? R:! 50 required for Yukawa unification, have an 
important impact on the SUSY spectrum and on the 
predictions for FCNCs. In particular, the branching 
ratio of the decay Bs — > fi'^l-t~ receives large tan/3- 
enhanced contributions from Higgs-mediated neu- 
tral currents that are proportional to (tan/3)^/M^ 
[IB, [i3| • With large tan /? and a large trilinear cou- 
pling At following from relation s fl] ), the stringent 
most recent experimental bound [l8| 

BR{B, M+/^")cxp < 5.8 X 10"^ (95% C.L.) (2) 

can only be met with quite heavy v4°, and i7+ 
Higgs bosons. 

Another important process in this respect is the 
tree-level decay t^v. Using the SM fit value 

for the CKM matrix element Vub [lO', '2^1 one obtains 
a SM prediction for the branching ratio^ 

BR(B+ ^ T+v)sM = (0.82 ± 0.11) x 10"" (3) 

that is quite low compared to the experimental value 

BR(B+ ^ T+i/)exp = (1-41 ± 0.43) X lO""* . (4) 



This prediction is obtained by normalizing the branching ra- 
tio to AMjj [2l[|. The value in eq. JSj agrees well with those 
reported in [gl. Il9l. [23|. 



In the MSSM with large tan/?, the dominant addi- 
tional contribution to this process comes from charged 
Higgs bosons and is found to interfere always destruc- 
tively [26| with the SM contribution, thus further 
reducing the theory prediction. Hence, similarly to 
Bs M^M J the _B+ — > T^v decay requires a heavy 
Higgs spectrum to be in agreement with the experi- 
mental data. However, given the large experimental 
uncertainty in ([¥]), the Bs ^J■^fJ■^ constraint turns 
out usually to be more stringent. 

Finally, a very important constraint is the inclu- 
sive decay B — > Xgj, which receives the dominant 
SUSY contributions from a chargino-stop loop and 
a top-charged Higgs loop.^ The chargino contribu- 
tion is tan /3-enhanced and, with the large negative 
trilinear parameters implied by relations ([T]), adds de- 
structively to the SM branching ratio. The charged 
Higgs contribution, on the other hand, adds construc- 
tively to the branching ratio, but is suppressed by 
the heavy Higgs masses required to be consistent with 
Bs — > fJ.'^ fi~ . Thus a near cancellation between the 
two contributions, which is necessary in view of the 
good agreement between the experimental determina- 
tion [27| 

BR(B ^ XsjUp 

= (3.55 ± 0.24t°:5-^ ± 0.03) x 10""* (5) 

and the SM prediction 28] 

BR(B -> A:^7)sm = (3.15 ± 0.23) x 10"'' , (6) 

is difficult to achieve. 

Note that the solution with the chargino contribu- 
tion so lar ge t hat the sign of the b ^ sj amplitude 
is reversed |29| is challenged in our framework by the 
experimental data on BR{B Xs(.+(.~) 
In fact, it would lead to a 3cr discrepancy between the 
prediction and the experimental figure for this branch- 
ing ratio. 

Leaving aside, for the moment, the possibility of a 
sign fiip in the 6 — > 57 amplitude (we will return to 
this issue in section Will , the above discussion implies 
that a tension between the prediction for B — > Xs'y 
and the bound on Bg — > /^^/i" should generally be 
expected in models with Yukawa unification, as a di- 
rect consequence of relations ([TJ and the large value 
of tan/3. By the nature of the argument, this tension 
should be completely independent of the mechanism 
(fiavour symmetries or other) embedded in the SUSY 
GUT to explain the light quark masses and mixings. 
In section HVl we will come back to this issue, show- 
ing that indeed this tension occurs generally in SUSY 



^ In the actual numerical analysis we include all the relevant 

contributions, in particular gluino-down squark loops as well. 

The latter are found to play a negligible role. 
^ This statement holds, barring non-negligible new physics 

contributions to the Wilson coefficients Cg'^Q (see Ref. [SQI), 

which is impossible in our case. 
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GUT models with Yukawa unification and quantifying 
the tension numericahy with a x^-procedure. 

Technically, the most immediate potential remedy 
to the above mentioned problem seems to be to lower 
tan/3. This in fact alleviates the pressure from Bg — > 
H'^ fj,~ , permitting in turn larger Higgs and smaller 
chargino contributions to _B ^ X^j and thereby mak- 
ing possible that those two contributions indeed can- 
cel to a large extent. 

Lowering tan /3 means breaking the unification of 
the top and bottom Yukawa couplings, so that the 
full Yukawa unification is relaxed to b — t unification, 
occurring, e.g., in SU(5). Such breaking oft — b unifi- 
cation is actually also a general possibility in SO(IO) 
SUSY GUTs once all the representations needed for a 
realistic GUT-breaking sector are taken into account. 
For example, the "minimal breaking scheme" intro- 
duced by Barr and Raby [s^l requires the presence 
of a 16'jLf spinor. In this framework, the MSSM Higgs 
doublet Hd can naturally be a mixture between a dou- 
blet contained in the same IOh representation as the 
doublet Hu and one doublet contained in this 16'^ 
spinor, since the two have the same quantum num- 
bers. One then has 0, i, Hi] 

- H{10h) , 
Hd = i7(10H)cos7 + i?(16'ff)sin7 . (7) 

Consequently, the Yukawa unification relation At — Xb 
is effectively broken to 

T^=cos7. (8) 
M 

At the EW scale, this relation leads to a value of 
tan /3 < 50 parametrically smaller than in the exact 
unification case, depending on the amount of mixing 
in the second of eqs. 0. 

We would like to emphasize that the two cases 
of SU(5) and SO(IO) with minimal breaking scheme 
mentioned above are just intended as examples. Our 
analysis will be completely general in SUSY GUTs 
with b — T unification. 

It should be stressed as well that, even without t — b 
unification, SUSY GUT models with 6 — t unification 
maintain in fact most of their predictivity, since the 
relation between the b and r Yukawa couplings re- 
quires the ISMH relations, eq. ([T]), to be satisfied 
in order to obtain a correct prediction for mi,. In 
addition, a crucial observation is that b ~ t unifi- 
cation requires tan/3 either close to unity (which is 
however excluded by the Higgs mass bound [1^) or 
0(50), because otherwise the predicted bottom quark 
mass is in general too large |36l . Although the 

case tan/3 = 0(50) can be significantly modified by 
the tan /3-enhanced threshold corrections to men- 
tioned above, 6— t unification is difficult to achieve for 
tan/3 < 35. Therefore the strategy to lower tan/3 is 
not a trivial one in our context, since b — r unification 
pushes by itself tan/3 to high values. 



With the above arguments, departure from third 
generation Yukawa unification and restriction to 6 — r 
unification seems to be a promising approach to retain 
the predictivity of GUT models, while at the same 
time possibly removing tensions in FCNC observables, 
thanks to tan (3 < 50. The rest of our Letter is thus 
an attempt to address the following two questions: 

• Is the tension between FCNC observables a gen- 
eral feature of GUT models with third gener- 
ation Yukawa unification and universal masses 
for sfermions and gauginos at the GUT scale; 

• Is this tension relieved when tan/? is (slightly) 
below 50, i.e., if one moves from exact Yukawa 
unification but retains b — r unification. 

These issues will be studied through a numerical pro- 
cedure to be described below. 



III. PROCEDURE 

We assume, at scales higher than the GUT scale 
Mc, the existence of a grand unified group entail- 
ing b ^ T unification. Beneath Mq the grand unified 
group is broken"* to the SM group Gsm = SU{3)c x 
SU{2)l X U{1)y and the MSSM running is performed 
down to the EW scale. As for the GUT scale bound- 
ary conditions to this running, we include a unified 
gauge coupling ac, allowing for a GUT scale thresh- 
old correction £3 to the strong coupling constant, the 
Yukawa couplings for up- and down-type fermions of 
the third generation At and Af,,^ and a right-handed 
neutrino mass M/j. At Mq we also assume a soft 
SUSY-breaking sector, consisting of a universal tri- 
linear coupling Aq, a universal sfermion mass mig, a 
universal gaugino mass mi/2, as well as non-universal 
Higgs mass parameters mn^ , rnn^ ■ 

We run all the parameters using one-loop RGEs for 
the soft sector and two-loop RGEs for the Yukawa 
and gauge couplings [s^. To take correctly into ac- 
count the effects of right-handed neutrinos present in 
SO(IO) and required for the see-saw mechanism, we 
include the contribution of a third-generation neu- 
trino Yukawa coupling (with initial condition A^ = 
At) in all RGEs between Mq and Mr [H, Ell- 
in our framework, there are thus no potentially large 
logarithmic GUT scale threshold corrections to either 
Yukawa unification or to Higgs splitting, which would 
be present if such contribution were neglected in the 
RGEs.^ The remaining GUT scale threshold correc- 



* For our purposes, this can be assumed to happen in one single 
step. 

^ On the Yukawa coupUngs of the Ughtest two generations we 

will comment later on in this section. 
^ In contrast with statements made in the literature and in 

accord with the results of [lOl | , we find that neutrino Yukawa 
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tions to the Yukawa couplings are expected to be small 

At the EW scale, we finally have the two additional 
free parameters tan f3 and fi. The total number of free 
parameters is then 13 and they are collected in table[Tl 



Sector 


# 


Parameters 


gauge 


3 


aa, Ma, 


SUSY 


5 


mi6, mi/2, ^0, mn,,, mn^ 


Yukawas 


2 


At, \b 


neutrino 


1 


Mr 


SUSY (EW scale) 


2 


tan/3, ^ 



TABLE I: Model parameters. Unless explicitly stated, 
they are intended at the GUT scale. 

After calculating the SUSY and Higgs spectra^ and 
the threshold corrections to third generation fermion 
masses |42| , we evaluate the flavour-changing observ- 
ables using the effective Lagrangian approach of [47| . 
Thereafter, in order to have a quantitative test of the 
model, we construct a function defined as 



(9) 



j-{ {(^t )cxp + {cFi )theo 

composed of the quantities given in tables HI] and HlI 



Observable 


Valuc((Tcxp) 


Observable 


Lower Bound 


Mw 


80.403(29) 


Mho 


114.4 


Mz 


91.1876(21) 




104 


lO^G^ 


1.16637(1) 


Mi 


95.7 


l/Qem 


137.036(0) 






as{Mz) 


0.1176(20) 






Mt 


170.9(1.8) 








4.20(7) 






Mr 


1.777(0) 







TABLE II: Flavour conserving observables [48l. |49| used in 
the fit. Dimensionful quantities are expressed in powers 
of GeV. 



In eq. (O Oi indicates the experimental value of the 
observables listed in tables [TTland lllll and fi\'ff\ the cor- 
responding theoretical prediction, which will be func- 
tion of the model parameters listed in table HI collec- 
tively indicated with d. The function is minimized 
upon variation of the model parameters, using the 
minimization algorithm MIGRAD, which is part of the 
CERNlib hbrary [53]. 

Some comments are in order on the determination 
of the errors. First, one can note that among the 



Observable 




Value(crcxp)(crthco) 


Ref. 


AM /AM^ 
lO"* BR(B ^ Xsi) 
10^ BR(B ^ XsC+r) 
10'' BR(B+ ^ r+i/) 
BR(B, ^ ^1+^.-) 


35.1(0.4)(3.6) 
3.55(26) (46) 
1.60(51)(40) 
1.41(43)(26) 
< 5.8 X 10"* 


[27, 50] 
[27] 

[51, 52] 
[23] 
[18] 



TABLE III: Flavour-changing observables used in the fit. 



The BR(B 
[1,6] GeV^ 



Xst^ t ) is intended in the range £ 



observables in table [TTl some have a negligible exper- 
imental error. In this case, we took as overall un- 
certainly 0.5% of the experimental value, which we 
consider a realistic estimate of the numerical error as- 
sociated with the calculations. Concerning the theo- 
retical errors on the flavour observables (table we 
note the following: the error on AMg/ AM^ takes into 
account the SM contribution, dominated by and 
the new physics contributions, dominated by scalar 
operators; the error on BR(i?+ t^v), after nor- 
malization by AMd pT|, is dominated by the lattice 
"bag" parameter Bd and the relevant CKM entries; 
the error on BR(i3 Xg^) is taken as twice® the 
total theoretical error associated with the SM calcu- 
lation [H]; finally the error on BR(B Xsi+i~) is 
taken as 25% of the experimental result, and is esti- 
mated from the spread of the theoretical predictions 
after variations of the scale of matching of the SUSY 
contributions. 

In evaluating the function, we also included the 
bounds reported in tables |TT] and IIIII These con- 
straints are in the form of suitably smoothened step 
functions, which are added to the x^-function of eq. 

If any of the constraints is violated, the step func- 
tions add a large positive number to the , while for 
respected constraints the returned value is zero, so 
that the is set back to its 'unbiased' definition ((9]). 

A step function was also included in order to enforce 
the desired sign for the 6 — > 57 calculated amplitude, 
thus permitting to systematically explore both cases 
of like sign or flipped sign with respect to the SM one. 
In the case of flipped sign, large SUSY contributions 
are necessary such that SUSY is not quite a correc- 
tion to the SM result, but rather the opposite. As a 
consequence, one would need a theoretical control on 
the SUSY part at least as good as that on the pure 
SM calculation. In the absence of this knowledge, the 
amplitude in the flipped-sign case is generally very 
sensitive to variations of the matching scale, and the 
associated theoretical error hard to control. In or- 
der to be able to estimate as reliably as possible the 



effects are not sufficient to explain the lar ge H iggs splitting 
required for successful EWSB to occur (see [id ], footnote 15). 
We calculate the Higgs VEVs and following \^ and use 
FeynHiggs [11,143,13,13 to obtain the light Higgs mass. 



This choice is quite conservative, considering that in our case 
variations of the calculated BR(B — > ^^7) upon variation of 
the SUSY matching scale in the huge range [0.1, 1] TeV are 
typically around 4%. However, we feel it is justified in the 
case of large cancellations among new physics contributions. 
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6 ^ S7 amplitude also in the flipped-sign case, we 
have taken advantage of the SusyBSG code [sj] , which 
is directly called by the fitting procedure. 

As already mentioned in section [Til the scenario 
with flipped 6 ^ 57 amplitude leads to ^ 9 
solely on account of the implied 3a discrepancy in 
B Xsi+i-. We calculated the BR{B Xst+f.') 
using the results of Ref. [HHl . We will address the sce- 
nario with flipped b ^ sj amplitude quantitatively in 
section IIVI 

An important observation is in order at this point, 
justifying why our analysis should be valid for any 
SUSY GUT model with b — r unification and universal 
soft terms as in table HI As already mentioned in sec- 
tionini any such model prefers the region of parameter 
space leading to ISMH, i.e., third generation sfermion 
masses much lighter than first and second generation 
ones, the latter being of 0(mi6). One should as well 
consider that, assuming hierarchical Yukawa matrices 
and large tan (3, it is sufficient to include only the 33- 
elements of Yukawa matrices in the RGEs, that is, 
take GUT-scale boundary conditions for the Yukawa 
couplings as 

Yu,d = diag(0, 0, Xt,b) ■ (10) 

This will have a negligible effect on the determination 
of the first and second generation sfermion masses, 
given their heaviness. This observation makes it pos- 
sible to separate the effects of specific Yukawa tex- 
tures, which depend on the flavour model one embeds 
into the SUSY GUT, from those genuinely due to the 
unification of Yukawa couplings. The adoption of this 
strategy brings us to the following remarks: 

• Given the approximation we adopt for the initial 
conditions on the Yukawas, we do not need to 
assume any particular flavour model. The low- 
energy input of the CKM matrix and of the fer- 
mion masses other than third generation ones, 
necessary for the calculation of many among the 
observables included in the fit, is then taken di- 
rectly from experiment^. 

• Due to ISMH, the lighter stop is always the 
lightest sfermion and in fact its tree-level mass 
can be very small. Therefore, we include the 
one-loop corrections to the light stop mass to en- 
sure our solutions are consistent with the lower 
bound in table Ull In practice, due to its light- 
ness, we calculate the stop pole mass with the 
same accuracy as the pole masses of the third 
generation fermions, the W, Z, and the Higgs 
bosons. 



^ Specifically, the CKM input is taken from the new physics 
independent CKM fit [Tj. 



• The validity of our approximation, eq. PU)) . was 

checked by performing the full analysis also with 
a SUSY GUT model with specific fiavour tex- 
tures, namely the DR model [l3| (with Yukawa 
unification relaxed as in eq. Our results 

were not significantly affected. 

IV. RESULTS 

In order to address the questions outlined in the 
introduction, we have explored in the mig vs. Xt/Xb 
plane the class of GUT models defined in section [TTTl 
Concretely, we fixed rriig to values > 4 TeV and Af / Xb 
to values > 1 (roughly equivalent to fixing tan (3 to 
values < 50) and minimized the function Q upon 
variation of the remaining model parameters. The 
minimum value provides then a quantitative test 
of the performance of the model. The results of our 
survey are reported in the four panels of Fig. [T] In 
particular, panels (a) to (c) report the contours as 
solid lines in the mig vs. tan (3 plane. As reference, 
also the values of Xt/Xb are reported on a right-hand 
vertical scale. Superimposed to the x'^ contours are: 
in panels (a) and (b), the deviations of respectively 
BR(i3 — > ^^7) and rub from the central values in ta- 
bles imillll in units of the total error; in panel (c) the 
mass contours of the lightest stop. Finally, panel (d) 
shows the contributions from BR{B — > ^^7), nib 
and all the rest (as three stacked contributions, rep- 
resented by solid lines) vs. mie in the special case 
of Ai/Ab = 1, corresponding to exact Yukawa unifica- 
tion. 

Various comments are in order on these plots. 

1. Panel (d) shows that, for any mig < 9 TeV, the 
X^ contribution from B Xgj alone is no less 
than roughly 4, corresponding to no less than 2a 
deviation from the result of eq. ([5]). Therefore, 
in the case of Yukawa unification, agreement 
among FCNCs is only achieved at the price of 
decoupling in the scalar sector. One should 
note in this respect the quite conservative choice 
of the B —> Xsj error, already mentioned in 
footnote [HI The apparent non-monotonic be- 
haviour of the B — > Xsj x^-pi'ofile is due to 
the fact that, for toiq < 6 TeV, the model pre- 
diction becomes so bad that the minimization 
algorithm prefers to sacrifice the prediction for 
rUb (which in fact gets much worse). 

2. For mi6 < 4.7 TeV, fits usually prefer the 
flipped-sign solution for the & — > 57 ampli- 



" Of course such test cannot be attached a statistically rigorous 
meaning, since, e.g., the x^-entries are not all independently 
measured observables. 

^ For similar findings in the context of Bayesian analyses of 
the CMSSM, see, e.g., Ref. 
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FIG. 1: Panels (a)-(c); ^ contours (solid lines) in the mie vs. tan/3 plane. Superimposed as dashed lines are the pulls 
for BR(_B — > y^s'i) (panel (a)) and for m;, (panel (b)) and the lightest stop mass contours (panel (c)). Panel (d): x'^ 
contributions vs. mis in the special case of exact Yukawa unification. All the plots assume a SM-like sign for the fo — > 57 
amplitude, except for panel (d), where also the total for the hipped-sign case is shown as a dot-dashed line. 



tude, discussed in section Hill However, at the 
quantitative level, this requirement (implying a 
quite light stop mass) turns out to be difficult 
to reconcile with that of successful predictions 
for the bottom mass and/or for EW observ- 
ables. The solutions we found have ^ 14 
as result of the contributions from EW observ- 
ables and rvih alone. On top of it one has to 
add the contributions from FCNCs, being 
> 9 solely on account of the 3tT discrepancy in 
BR(B Xs£+i-). The total x^ for the fits 
with flipped-sign b ^ sj amplitude is reported 
in figure [11(d) as a dot-dashed hue. 



3. From panels (a)-(c) one can note a region of 
successful fits for mig > 7 TeV and 46 < 
tan/3 < 48, corresponding to a moderate break- 
ing of t — b unification, since it corresponds to 
0.2 > (At/Ab - 1) > 0.1. 

4. The interesting region, characterized in point|31, 
emerges as a compromise between B — > Xs"f 
and rrif,, pushing tan /3 to respectively lower and 
larger values. In this region, the lightest stop 
mass is below order 1 TeV, but not less than 
roughly 800 GeV. 

5. If mig is not too large, the interesting region 
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90.73 
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171.3 


0.2 


mb(m6) 




4.20 


4.28 


1.1 






1.777 


1.77 


0.4 


10* BR(B - 




3.55 


2.72 


1.6 


10^ BR(B 


1.60 


1.62 


0.0 


AMs/AMa 




35.05 


32.4 


0.7 


10* BR(B+ 


^ T+U) 


1.41 


0.726 


1.4 


10** BR(Bs 




< 5.8 


3.35 








total X'''- 


8.78 





Grs 


Mass predictions 


mi6 


7000 


MhO 


121.5 




1369 


Mjjo 


585 


Ml/2 


143 


Ma 


586 




-14301 


Mh+ 


599 


tan /3 


4D.i 


■mil 


783 


1/«G 


0/1 7 




1 TOO 

1 (28 


iVlQ 1 lU 


o.u / 








-4.91 




2378 


(rriu l'm^a'\^ 


1.616 




3297 


{mHa/mief 


1.638 




58.8 




8.27 




117.0 




0.608 




117.0 


Ad 


0.515 




470 



TABLE IV: Example of successful fit in the region with b — t unification. Dimensionful quantities are expressed in 
powers of GeV. Higgs, lightest stop and gluino masses are pole masses, while the rest are running masses evaluated at 
Mz. 



Observable 




Exp. 


Fit 


Pull 


Mw 




80.403 


80.32 


0.2 


Mz 




91.1876 


90.63 


1.2 


10^ G„ 




1.16637 


1.162 


0.7 


I/Cfcm 




137.036 


136.4 


1.0 


as{Mz) 




0.1176 


0.1144 


1.5 


Mt 




170.9 


171.5 


0.3 


mb{mb) 




4.2 


4.41 


2.9 


Mr 




1.777 


1.77 


0.7 


10* BR(B ~* Xsj) 


3.55 


3.69 


0.3 


10^ BR{B - 




1.60 


4.41 


4.3 


AMs/AMd 




35.05 


32.5 


0.7 


10* BR(B+ 




1.41 


0.561 


1.7 


10* BR(Bs 




< 5.8 


5.0 








total x''- 


36.5 



Input parameters 


Mass predictions 




4000 


M„o 


114.5 




249 


Mjjo 


509 


Ml/2 


149 


Ma 


510 




-7989 


M„+ 


519 


tan /3 


50.2 




425 


1/qg 


24.7 




823 


Mg/IQ^^ 


2.24 




680 


£3/% 


-3.78 




759 




1.643 




1402 




1.908 


^1 


59.3 


Mfl/lO^^ 


3.30 


X2 


109.3 




0.643 




108.9 


Ad 


0.643 


Ms 


461 



TABLE V: Representative fit obtained in the region with t — b — r unification and flipped b - 
as in table HVl 



S7 amplitude. Conventions 



is clearly distinguished from the corresponding 
case with exact Yukawa unification, as far as 
the fit quality is concerned. By looking at pan- 
els (a)-(c), one can in fact recognize that the 
gradient of variation is close to vertical for, 
say, mi6 = 7 TeV, when increasing tan /3 from 
around 48. On the other hand, increasing niie 
makes the case of breaking t — b unification more 
and more indistinguishable, for the fit perfor- 
mance, from the decoupling regime of point 1. 

6. An example of a fit in the interesting region is 
reported in table lIVl and shown in panels (a)-(c) 
as a black square. Note that the prediction for 
BR(i3 Xs'y) still tends to be on the lower side 
of the experimental range from eq. ([5]), with a 
central value around 2.9 x 10""*. As a compari- 
son, a representative fit in the region with exact 
Yukawa unification, featuring a fiipped b sj 
amplitude, is reported in table El 

We note, point 1. above implies that, for any 
SUSY GUT with Yukawa unification, compatibility 



with FCNC observables and with the now precisely 
known value of nib selects the "partially decoupled 
solution" , advocated in Tobe and Wells [H^l , as the 
only phenomenologically viable. The low-mie solu- 
tion originally found in p^ . [l3l | is disfavoured when 
combining all the most recent data. This conclusion 
has been here quantitatively assessed with a pro- 
cedure.^^ 

Points [31 to [HI illustrate instead that compatibility 
among all the considered observables at values of mie 
of order 7 TeV can be recovered at the price of relaxing 
t—b—T unification to just b—r unification and without 
modifying universalities in the soft terms at the GUT 
scale. 



Note that, in our approach, exact Yukawa unification can 
be enforced, so that the lower bound on mig emerges trans- 
parently as a tension among observables. Instead, in, e.g., 
[58j . it is low-energy observables (like rUb) to be fixed and a 
large value for mig is needed for Yukawa unification to occur 
within a given tolerance. 
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We conclude this section with a few additional re- 
marks. A first interesting issue is whether the gen- 
eral tension involving FCNCs and nib studied in this 
Letter may be relieved if, instead or in addition to 
lowering tan /?, one allows for a complex phase in the 
trilinear coupling At . The latter in our analysis is real, 
since we take real Aq. A complex phase 0t in At would 
induce a cos suppression factor in the leading char- 
gino correction to BR(i? — + ^^^7). In the presence of 
a complex At, however, the tan /3-enhanced chargino 
contribution to the bottom mass would also become 
complex, so that a (chiral) redefinition of the down- 
quark fields would be necessary in order to end up 
with real and positive masses. This may impact non- 
negligibly the overall size of SUSY corrections to mt 
Il4| as well as the pattern of CP violation in B- 
physics observables. In addition, a complex phase in 
At is also constrained by the electric dipole moments 
of the electron and neutron. Addressing these issues 
quantitatively goes beyond the scope of the present 
Letter, where we confine ourselves to real GUT-scale 
soft terms. 

As a second remark, we note that, in our analysis of 
the parameter space, we restricted ourselves to posi- 
tive values of the fi parameter. This is because, for 
negative fi, the ISMH solution is lost @, thus lead- 
ing to much heavier third generation sfermions for a 
given value of mig (typically > 2.5 TeV), which is 
unattractive for naturalness reasons. In addition, this 
would lead to a negative contribution to the anoma- 
lous magnetic moment of the muon a^, whereas a 
sizeable positive contribution is currently favoured by 
experiment. 

In fact, in the region preferred by the function, 
the contributions to are much smaller than would 
be needed to explain the E821 result [H,!!^!- This is a 
well-known fact [1, |Tl| and we decided not to include 
this observable in our function to obtain results 
unbiased by this issue, which still needs to be settled. 

We also chose not to include the constraint from the 
'standard' neutralino relic density calculation. The 
latter yields a posteriori a too large relic density in the 
interesting region of parameter space for well-known 
reasons, i.e., the heaviness of sfermions and that of 
the pseudoscalar Higgs, making the A-funnel region 
unattainable. This constraint can however be circum- 
vented in many ways, such as a late-time entropy in- 
jection (see, e.g., [61|).^'^ 

For the above reasons, we preferred to restrict the 
set of observables considered in the analysis to those 
listed in tableslTTland lllll on which consensus is broad. 



The requirement that SO{10) SUSY GUTs with Yukawa uni- 
fication do account in full for the standard CDM abundance 
has been recently reconsidered in j63 |. 



V. CONCLUSIONS 

In this Letter we have studied the viability of the 
hypothesis of t — & — r Yukawa unification in SUSY 
GUTs under the assumption that soft-breaking terms 
for sfermions and gauginos be universal at the GUT 
scale. We found that this hypothesis is challenged, 
unless the squark spectrum is pushed well above 1 
TeV. Our conclusion is assessed through a global fit 
including EW observables as well as quark FCNC pro- 
cesses. The origin of the difficulty is mostly in the spe- 
cific parameter region chosen by Yukawa unification, 
which guarantees the correct value for the bottom 
mass and implies ISMH. In this region, it is impossi- 
ble to accommodate simultaneously the experimental 
value for BR(i3 —>■ Xs^) and the severe upper bound 
on BR(i?s — >• fi^ fjL'). This statement holds under the 
prior assumption that the sign of the b sj ampli- 
tude be the same as in the SM. For mig < 4.7 TeV, 
fits prefer the flipped-sign solution for the 6 — > 57 
amplitude. In this instance it is however difficult to 
achieve agreement on the 6 — > 57 prediction and, si- 
multaneously, on those of EW observables and/or of 
the bottom mass. In addition, this possibility entails 
in our case a 3<t discrepancy in B Xs£~^£~ data. 

We have shown that our above conclusions hold 
irrespectively of the flavour model embedded in the 
SUSY GUT, since, due to ISMH, first and second 
generation squark masses are much heavier than third 
generation ones. 

We have also addressed the possibility of relaxing 
Yukawa unification to 6 — t unification, which in fact 
allows to maintain most of the predictivity of SUSY 
GUTs. In this case, the tension underlined above is in 
fact largely relieved. The fit still prefers large values 
of 46 < tan /? < 48, as a compromise between FCNCs 
and 7716, pushing tan /3 to respectively lower and larger 
values. The range for tan P corresponds to a moderate 
breaking of i — 6 Yukawa unification, in the interval 
At/Afc-1 e [0.1,0.2]. 

In the interesting region, we find the lightest stop 
mass > 800 GeV, a light gluino around 400 GeV 
and lightest Higgs, neutralino and chargino close to 
the lower bounds. This spectrum implies BK{Bs 
fi'^fi^) in the range 2 to 4x10"^ and BR{B -> X^'j), 
on the lower side of the "acceptable" range, « 2.9 x 
10~*. We stress that the requirements of 6 — r unifi- 
cation and the FCNC constraints are enough to make 
the above figures, exemplified in table IIV) basically 
a firm prediction within the interesting region, hence 
easily falsifiable once the LHC turns on. 
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